Twirling and Whirling: Viscous Dynamics of Rotating Elastica 
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Motivated by diverse phenomena in cellular biophysics, including bacterial flagellar motion and 
DNA transcription and replication, we study the overdamped nonlinear dynamics of a rotationally 
forced filament with twist and bend elasticity. Competition between twist injection, twist diffusion, 
and writhing instabilities is described by a novel pair of coupled PDEs for twist and bend evolution. 
Analytical and numerical methods elucidate the twist/bend coupling and reveal two dynamical 
regimes separated by a Hopf bifurcation: (i) diffusion-dominated axial rotation, or twirling, and 
(ii) steady-state crankshafting motion, or whirling. The consequences of these phenomena for self- 
propulsion are investigated, and experimental tests proposed. 

PACS numbers: 87.16.-b, 46.70.Hg, 47.15. Gf, 05.45.-a 



Dynamics and stability of rotationally forced elastic fil- 
aments arise in several important biological settings in- 
volving bend and twist elasticity at low Reynolds num- 
ber. In the context of DNA replication, when two daugh- 
ter strands are produced from a duplex, it was noted [Q 
long ago that energy dissipation for rotations about the 
filament axis is so much smaller than that for transverse 
motions that axial "speedometer-cable" motions are fa- 
vored, and are energetically and topologically feasible. 
During DNA transcription, in which a polymerase pro- 
tein moves down the double-stranded filament, progres- 
sive unwinding of the helix can lead to an accumulation 
of local twist that may induce "writhing" instabilities 
of the filament ||. Energetic and dynamical aspects of 
these processes are of great current interest |M. 

At the cellular level, bacteria are propelled through 
fluids by helical flagella turned by rotary motors in the 
cell wall ||. Recent studies ||] have revealed the details 
of two competing crystal structures assumed by flagellin, 
the protein building block of flagella, corresponding to 
helices of opposite chirality. Both local and distributed 
torques can change the conformation of flagella; during 
swimming these motors episodically reverse direction , 
and the resultant torques can induce transformations be- 
tween these states ||, while uniform flow past a pinned 
flagellum may induce such chirality inversions 

To elucidate fundamental processes common to these 
systems, we consider here the model problem shown in 
Fig. H a slender elastic filament in a fluid of viscosity r\, 
rotated at one end at frequency ujq with the other free. 
We study competition between three processes: twist in- 
jection at the rotated end, twist diffusion, and writhing. 
Analytical and numerical methods reveal two dynamical 
regimes of motion: twirling, in which the straight but 
twisted rod rotates about its centerline, and whirling, in 
which the centerline of the rod writhes and crankshafts 
around the rotation axis in a steady state. 

This work is a natural outgrowth of recent studies of 



forced elastica in the plane |lC||Ll|j , and dynamic twist- 
bend coupling p^-p^|. The balance considered between 
elastic and viscous stresses complements that between 
elasticity and inertia in the inviscid limit (as in whirling 
shafts p5|Jl^ ]), where twist waves propagate ]l5|Jl7| . 

An clastic filament is characterized by its radius a, con- 
tour length L, bending modulus A, and twist modulus C. 
The total elastic energy cost £ for curvature k and twist 
density f2 is an integral over arclength s parameterizing 
the position r(s,t) of the filament centerline EBI, 
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where the Lagrange multiplier A enforces local inextensi- 
bility, (rt) s ■ r s — 0. Thus arise two dimensionless ratios: 
r = C/A, and the aspect ratio a/L. At zero Reynolds 
number, elastic forces per length f = —5£/5r balance 
the viscous drag from slender-body hydrodynamics fll9|: 



C||ttT t + a(I-tt)-r 4 = f, 



(2) 



where t = r s is the unit tangent, and the transverse 
and longitudinal drag coefficients are (± ~ 2£|j ~ 
4irrj/[hi(L/2a) + c], with c a constant of order unity Jl9| . 
Likewise, the axial elastic torque per unit length m = 
Cfl s Jl4|,[l5| balances the local rotational drag: m — C, r uj, 
where uj(s,t) is the local rotational velocity about t and 




FIG. 1. An elastic filament, rotated about z at the left end, 
surrounded by a fluid of viscosity n. 
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Cr 47rr/a 2 @. We also define e 2 = Cr/(G-£ 2 ), so 
e ~ (a/L), apart from logarithmic corrections. 

The dynamics are closed by a geometric constraint, 
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which shows how twist changes due to differential rates of 
angular rotation, stretching (e.g., extension of a straight, 
twisted rod decreases f2), and out-of-plane bending mo- 
tions (writhing) Q-0| . The constraint (pf) is a conserva- 
tion law for twist density f2, with twist current —to, and 
with the stretching and writhing terms acting as sinks or 
sources. The velocities u> and r t enter (|J) through their 
space derivatives, since rigid motions cannot change Q. 
With the local torque balance C r uj = CQ S and assuming 
inextensibility we obtain 
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The second term of (|J) is nonzero |20|] when the filament 
is both out of elastic equilibrium (f ^ 0) and either non- 
planar (with torsion t ^ 0) or twisted (fi ^ 0), and then 
acts as a sink or source of twist. 

For boundary conditions, we assume the forced end 
(r(0) = 0) of the rod is clamped (r s (0) = z) and the free 
end experiences no force or torque (t ss (L) — r sss (L) = 
0; fi(L) = 0). Local torque balance sets f2 s (0) = C, r uJo/C . 

Before solving these PDE's numerically, we use dimen- 
sional analysis to understand the main features of the 
motion. We focus on small-amplitude bend and twist de- 
formations of a straight filament (thus ignoring Of 
the seven parameters (A, C, L, a, Cr, Cu w o), four remain 
after introducing Y and e and noting that a only appears 
through ( r , so suitable rescalings of length, time, and f2 
leave only one control parameter. This can be chosen 
proportional to the rotation frequency ujq. 

For low turning rates loq, the filament remains straight, 
with twist diffusion and injection balancing. The twist 
density at the clamped end follows from a balance of 
viscous and elastic twisting torques, Cr<-<JoL ~ CT2, or 
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Instability occurs when the the twist torque C'£l is com- 
parable to the filament buckling torque A/L (T3|. At this 
point, the balance of viscous and twist torquesflq) implies 
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where the second and third result follow from the rela- 
tions A = (it/4)a 4 E between the bending modulus and 
the Young's modulus E |15j and A = ksTL p , with L p 
the persistence length. Interestingly, lu c is independent 
of the twist modulus C, and since the twist density scales 
with the drag, lu c varies inversely with £ r . 

The result ([|) is central, for naive dimensional anal- 
ysis predicts ui c ~ E/r). With E ~ 10 7 dynes/cm 2 and 



T) ~ 0.01 Poise as for rubber in water, co c would be enor- 
mous if not for the prefactor (a/L) 2 . Since a/L is rea- 
sonably 10~ 3 , we find ui c ~ 10 3 s _1 , similar to flagella 
rotation rates p| and achievable in the laboratory. The 
rightmost form in (Q) shows readily the frequency scales 
for systems of varying length and stiffness. Consider the 
elastica DNA (L p ~ 5 x 10~ 6 cm, a ~ 10~ 7 cm), mi- 
crotubules (L p ~ 0.5 cm, a ~ 10~ 6 cm), and bacterial 
filaments (L p ~ 25 cm, a ~ 3 x 10~ 5 cm) 0. The fre- 
quencies kBT/( r L p are then 8 x 10 6 s _1 , 0.8 s _1 , and 
1.3 x 10~ 5 s _1 , respectively. Thus, for the instability to 
appear at, say, 10 3 s _1 requires a minimum ratio L/L p of 
90, 0.03, and 10" 4 . A strand of DNA with L/L p ~ 10 2 is 
clearly not straight in isolation, so this instability would 
be hard to realize in DNA, but the stiffer examples of mi- 
crotubules and bacterial filaments are indeed candidates. 

Linearizing (^) and (|J) about a straight filament along 
z, with r « sz + r± , we see that twist diffuses with diffu- 
sion constant C/( r , while the backbone obeys a "hyper- 
diffusion" equation Ci r _Lt = — ^r^ ss;5S + Cz x (Orj_ ss ) s . 
For crankshafting motions, we set r± t = \ z x r _i_- Thus 
we find two characteristic lengths [|l0| , 

t x (x) = (A/C±x) lfi and = (C/C,^ ) 1/2 . (7) 

These are analogous to the penetration depth in the fa- 
miliar theory of oscillations in a viscous fluid |22| . The 
primary instability is given by £ r (uj) ~ L. 

The crankshafting frequency \ f° r whirling can be es- 
timated by assuming that the transverse drag, C±xl r J-|> 
is roughly equal to the elastic force per length, A|rj_4 S | ~ 
A\r ± \/L 4 . Thus X ~ A/C, X L 4 ~ (a/L) 2 uj c (C ps A 
for typical materials p5|), and £±(x) ~ £ r (u c ) ~ L at 
the transition. The whirling rod does not undergo sim- 
ple rigid body rotation; the speedometer-cable rotational 
motion is faster than the backbone crankshafting motion 
by a factor of (L/a) 2 . This steady-state shape is pos- 
sible because diffusion can homogenize the twist as fast 
as backbone motion relieves it. We describe this process 
quantitatively by integrating (|J) along the rod for inex- 
tensible, steady-state (fij = 0) crankshafting. The differ- 
ence Aoj in rotational velocities about the local tangents 
at s = L and s = is 



Aw=x[l-*-t(L)]. 



(8) 



Equivalently, Aoj is the injected twist current minus the 
twist current leaving the free end. Thus, writhing acts as 
a twist sink in steady-state crankshafting motion when 
the rod's free end is not aligned with the z-axis. 

When the twist diffusion time, ( r L 2 /C, is longer than 
the bending time x _1 , buckling can relieve twist faster 
than it is replenished by diffusion, and steady-state 
crankshafting would likely be unstable. One possible new 
behavior would consist of repeated sequences of transient 
whirling followed by quiescence as twist builds up anew. 
Scaling arguments yield a critical frequency of luq ~ E/r), 
a factor of e~ 2 higher than the rate at onset of the first 
instability, and thus unreachable for typical materials. 



2 



The bend relaxation time suggests the rescaling, 

?=(A/ax 4 ) . (9) 

A natural pair of further rescalings of (Q) is s — s/L, and 
f2 = £IL. If we parameterize the filament centerline (Fig. 
|l|) asr(s,i) = (X(s, t),Y(s, t),s — 5(s,t)), introduce the 
complex transverse displacement £ = (X + iY)/L, and 
expand the dynamics up to third order in £ (immediately 
dropping the tildes), we obtain 

6 = -U, - (A&), + *r p - s s ) + 5 SS $ S )} S 
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where the inextensibility constraint (1 — £ii /£_[_) (r ss • f) + 
(r s • f s ) ~ ^(r ss • f) + ( r s • f«) = is expanded to set A as 
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This constraint also fixes i5 s ~ §|£ s | ■ As anticipated, 
apart from material properties T and e, the coupled 
twist/bend dynamics are governed by a single control pa- 
rameter loq, which appears only in the (rescaled) bound- 
ary condition, fi a (0) = e 2 Lu /T = a. The PDEs @ are 
like those of excitable media Q, with a separation of 
time scales derived from the aspect ratio e. For the usual 
case £Cl, twist is the fast variable and bends are slow. 

Linearized Dynamics: The twist profile in the straight 
filament (£ = 0) satisfies fit = (T/e 2 )fl ss . After tran- 
sients die out, the steady-state profile is linear in s ||, 
f2 = (Cr^ofC) (s — L) . Using this in the linearized fil- 
ament evolution and taking £(s,t) = £(s) exp(i^t), a 
rigidly rotating, neutrally stable shape, we obtain 



x£ = *£4. + ra[(a -!)£„]« 



(11) 



Numerical solution |24[] of (|TT| ) yields a critical value 
a c — 8.9/r (confirming dimensional analysis of Eq. ([)])), 
below which the rod is straight and executes only axial 
rotation ("twirling"), and above which the rod buckles 
and rotates ( "whirls" ) at a frequency which for a ~ a c is 
X ~ 2.32r<x This motion is the dynamical equivalent of 
the static writhing instability of a twisted elastica p5[ . 
Inserting all numerical factors, the critical frequency and 
rotation rate (at onset) are 
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Weakly Nonlinear Theory: Numerical solution of (10) 
with a pseudospectral method 0,^5) shows that there is 
indeed a steady state beyond the bifurcation. As a — a c 
increases, the shape becomes more helical. The free end 
of the whirling filament experiences more drag than 




FIG. 2. (a) Stroboscopic montage of the "whirling" fila- 
ment, viewed from along the z-axis, as it rotates clockwise; 
(b) side view at two times, drawn to a different scale. 

points closer to the driven end and thus lags behind 
(Fig. ||). Since fl depends quadratically on the back- 
bone shape (see Eq. |3|), near w c , where £ is small, the 
twist density remains nearly linear in s. Numerical stud- 
ies show that the free end traces out a circle with radius 
R ~ (luq — u^) 1 / 2 : a supercritical Hopf bifurcation p3j. 
This can be understood from (|8[) and dimensional argu- 
ments for the displacement; [l — z ■ t(L)l ~ (R/L) 2 , so 
w(0) — u>(L) ~ x{R/L) 2 ■ The twist that can be relieved 
by diffusion is limited, so U)(V) w oj c . From the linear 
dynamics, % ~ w o leading to R ~ (cjq — ^c) 1 ^ 2 - 
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FIG. 3. (a) Amplitude R (with a square root fit) and (b) 
crankshafting frequency of filament tip motion as a function of 
frequency offset from primary instability, (c) and (d); filament 
shapes for (u>o — uj c )/lo c — 0.27 and 3.52, with F = 1. 
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FIG. 4. Propulsive force generated by steady-state whirling 
motions, as a function of driving frequency. 

Swimming: Chirality of the whirling filament breaks 
time-reversal invariance of the motion, thereby allow- 
ing |26| a net propulsive force F p along z to be gen- 
erated. The elastic propulsive force density is a total 
derivative, so the total force is expressible in terms of 
the filament properties at its ends. For the clamped/free 
boundary conditions used here, F p — — ^(0) -z — A(0) = 
r ss (0) ■ r ss (0) — A(0). As shown in Fig. f|, F p rises lin- 
early from zero near the bifurcation, as it is quadratic 
in the transverse displacement, which in turn has the 
supercritical form shown in Fig. |^(a). While we know 
of no organism that utilizes this precise mechanism for 
self-propulsion, there is evidence in certain laboratory ex- 
periments pTj] for self-propulsion associated with twist- 
induced whirling in growing bacterial macrofibers con- 
strained at one end. Experiments are underway to ex- 
plore further this possible connection. 

The possibility of observing instabilities driven by twist 
accumulation along an elastica hinges upon a balance of 
material properties, fluid viscosity, and adequate forc- 
ing. Flow- and rotation-induced bacterial flagellar con- 
formational transitions provide proof-of-principle 
that this balance can be achieved in vivo. Like flagella, 
fibers of B. subtilis cells have adequate material proper- 
ties (e.g. Young's modulus and aspect ratio to dis- 
play instabilities like those described here. More complex 
phenomena are associated with instabilities of rotating 
helical flagella, as described elsewhere p8[ . 
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